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ABSTRACT 



We detect anisotropy in the cosmic microwave background (CMB) 
at degree angular scales and confirm a previous detection reported by 
Wollack et al. (1993). The root-mean-squared amplitude of the fluctuations 
is 44^7^ yuK. This may be expressed as the square root of the angular 
power spectrum in a band of multipoles between l^jf = 69^22- We find 



STi = ^l{2l + l) < |aj"|2 > /47r = A2tl^fiK. The measured spectral index of 
the fluctuations is consistent with zero, the value expected for the CMB. The 
spectral index corresponding to Galactic free-free emission, the most likely 
foreground contaminant, is rejected at approximately 3a. 

The analysis is based on three independent data sets. The first, taken in 
1993, spans the 26 - 36 GHz frequency range with three frequency bands; 
the second was taken with the same radiometer as the first but during an 
independent observing campaign in 1994; and the third, also take in 1994, spans 
the 36-46 GHz range in three bands. For each telescope position and radiometer 
channel, the drifts in the instrument offset are < 4 /iK/day over a period of one 
month. The dependence of the inferred anisotropy on the calibration and data 
editing is addressed. 

Subject headings: cosmic microwave background — cosmology: observations 
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1. Introduction 

In Wollack et al. (1993), we reported on a measurement of the anisotropy in the CMB 
made from Saskatoon, SK (SK93). The results were most parsimoniously described as a 
detection of CMB anisotropy. However, an anisotropy in Galactic free-free emission could 
have produced similar data with ^ 10% probability. We have since returned to Saskatoon 
and observed the same fields with an an extended frequency baseline that now spans 26 to 
46 GHz. The new data, SK94, confirm the 1993 results and strongly constrain the spectral 
index of the fluctuations. Reviews of the state of CMB anisotropy measurements are given 
in Bond (1994), Readhead & Lawrence (1992), and White et al. (1994). 



2. The Instrument and Observing Strategy 

Observations were made with two independent radiometers, one in Ka band and one in 
Q band, that sequentially shared the same telescope. The details of the instrument may be 
found in Wollack (1994a, 1994b), Wollack et al. (1993, hereafter Wl) and Page et al. (1994). 
Each radiometer detects the total power in three frequency bands and two polarizations^. 
The radiation is amplified by HEMTs (Pospieszalski et al. 1988; Pospieszalski 1992 ) 
cooled to 15K. The passbands for K^ band are 26-29, 29-32, 32-36 GHz, and for Q band 
are 36-39.5, 39.5-43, 43-46 GHz. The six channels in each radiometer view the sky in a 
single approximately frequency- independent beam formed by a cold corrugated horn feeding 
an ambient temperature off-axis parabola. The FWHM beam in Ka band is 1.42 ± 0.02° 
and in Q band is 1.04 ± 0.02°. After the parabola, the beam reflects off of a 90 cmx 150 



^A cracked waveguide joint vitiated the data from the vertical polarization in the Q 
radiometer; it is not used. 
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cm chopping flat that scans the beam back and forth about a vertical axis, following an 
approximately triangular waveform at 3.90625 Hz. A large stationary ground screen shields 
the radiometer from terrestrial and solar emission. Other than the new radiometer and 
the modulation scheme described below, the principle improvements to the instrument 
described in Wl are 1) an enlarged ground screen; in our data, contaminating signals 
from the ground and sun are less than 1 /iK. 2) thermal stabilization of all radiometer 
components to better than 0.1 K; and 3) increased vibrational isolation of the radiometers. 

The chopping fiat sweeps the beam in azimuth while the radiometer channels are 
synchronously sampled at 250 Hz. The the angle between the sweep endpoints subtends 
7.01° and 7.36° for the Ka and Q systems respectively; thus, the 3dB antenna beam 
contours of the two systems match at the extrema of their sweeps. The center of the sweep 
pattern is wobbled every 20 seconds from 4.42° west of the north celestial pole to 4.42° east. 

The effective beam pattern is synthesized in software by appropriately weighting 
each sample during a chopper sweep. This allows us to probe multiple angular scales, 
optimize the spatial frequency coverage to mitigate atmospheric contamination, examine 
the quadrature demodulation, and match the Ka and Q synthesized profiles even though 
the FWHMs are different. The full analysis based on this scheme will be discussed in a 
forthcoming paper; here, we use the method to synthesize the SK93 pattern. The beam 
pattern for one chopper sweep is shown in Figure 1. The data are normalized so that 
a 1 /iK temperature change in the positive lobe of a beam averaged over 15° in right 
ascension (1 bin) gives a signal of 1 /iK. A good approximation to the window function is 
Wl = [1.5 — 2Pi{cos[9t]) + .5Pi(cos[2^t])]e"^'+^^'^ where Pi are the Legendre polynomials, 
Ot = 2.57° and a = 0.64°. 

The instrument is calibrated using Cassiopeia-A. From a fit of the 8-250 GHz data in 
Baars (1979) and Metzger (1986), we find S^ = (2070 ± i62)z/(-0-695±0.029) j^^ (^ ^^ q^Lz, 
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epoch 1994) corresponding to a 9.57 ± 1.2 niK signal with a 1.04° beam at 38.3 GHz. The 
laboratory calibrations agree with this to within 20%. Because Cas-A is observed in the 
same thermal and radiometric environment in which the CMB is observed, we consider it 
the definitive calibration. Cas-A is also used to map the beam; the values are within 0.02° 
of those found with a microwave source. The absolute pointing error, as determined from 
Cas-A, is 0.06°. 



3. Data Reduction and Analysis 

In essence, the data reduction consists of editing as described below, binning the 
results on the sky, and subtracting the mean of the binned data. The results are shown 
in Figure 2. Although the constant part of the demodulated signal, or offset, depends on 
the east/west position and on the frequency and polarization, it is relatively independent 
of time. A typical drift is 3 ± 3 yuK/day in both radiometers. Variations in the chopping 
flat temperature will effect the offset because of the dependence of the emissivity on the 
incident angle. The offset is fltted to the temperature of the chopper and a correlation is 
found; however, removing the small variation has a negligible effect on the final data set. 

There are three components to the data editing procedure. 1) The data from the 
single-difference demodulation Q, which is sensitive to horizontal spatial gradients in the 
sky, are used to monitor the atmosphere. If the mean deviation of 8s averages of data is 
greater than some value, rj, then the 15 minute section containing those data is cut. In 
Table 1, we have converted rj into the standard deviation of the variation in the horizontal 



^ "single- difference" refers to the demodulation where the data taken with the chopper 
pointed west are subtracted from those where the chopper is pointing east. 
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thermal gradient of the atmosphere, (. This value is relatively independent of the observing 
strategy. In addition, for a 15 minute section to be included, both abutting 15 minute 
sections must also be good. 2) Next, the standard deviation of the 64 datum from each 
chopper cycle is computed. If one or more points deviates by 3.5o", then the entire chopper 
cycle is deleted. This cuts an additional 2.6% of the good data. 3) Finally, we require that 
no 20s average of data in a right-ascension bin deviate by more than 3.5o" of all the data in 
the bin. This cut rejects ~ 1.5% of the data that passed the atmosphere cut and the final 
results are robust to variations in this cut level. 

The data quality is monitored with a quadrature demodulation (data taken during the 
clockwise chopper scan minus data taken during the counterclockwise scan). This signal is 
consistent with zero. 

To quantify the anisotropy signal, we use the Bayesian method described in Wl to 
find the root-mean-squared amplitude. Arms and index, f3, of the fluctuations. In the full 
implementation, the new data are analyzed along with the SK93 data using a 720 x 720 
correlation matrix. A maximum of the likelihood, L{Arms, P), is found when A^ms and 
f3 scale the fixed spatial correlations inherent in the observing pattern so that they best 
describe the patterns in the measured data. The data taken in the east and west are 
treated independently and the atmosphere and HEMT induced correlations are accounted 
for (Dodelson & Kosowsky 1994a). The resultant A^ms depends on the observing strategy. 
We convert it into an estimate of square root of the angular power spectrum in a band of 
multipoles between / = 47 and / = 98 (/g// = 69^22 ) following Bond (1994) and Peebles 
(1994). The results are given in Table 1. 
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4. Discussion 

Table 1 shows that there is a strong detection of anisotropy with (3^0, independent 
of atmosphere cut and data sub-set. If the Ka93 data are subtracted from the Ka94 data, 
the result is consistent with zero, as expected if both see the same signal. In other words, 
the Ka94 results confirm the Ka93 results. 

There are two dominant sources of systematic error in this measurement; inaccuracy 
in our knowledge of the fiux from Cas-A (13%) and an inaccuracy in our measurement 
of Cas-A (Q: 6%, Ka: 3%). Small errors are introduced by inaccuracies in the beam 
FWHM (2%), in the measurements of the passbands (2%), and in the phases of the 
recorded signals with respect to the optical axis (< 1%). The combination of all these 
errors leads to a ±15% uncertainty in our temperature scale. For the best estimate of 
celestial fluctuations we give the result with the most statistical weight (44^6^ /uK) with 
the ±15% added in quadrature. Arms = 44^7^ /iK. Converting to a band-power^, we find 



6Ti = Jl{2l ± 1) < laf^P > /47r = 42t^>K at kff = 69 



12,, T^ „+ ; _ «n+29 
-22- 



In Table 1, one notes that the addition of new data did not improve the statistical 
error bar by the square root of the observing time. In part, this is because the error on 
Arms contains the sample variance, which is a function of the area of the sky observed. The 
statistical error bar for the Ka93 data is ~ 10.5 /iK (the expectation value, opposed to the 
maximum of the likelihood, is roughly half way between the 16% and 84% integrals). If the 
sample variance, which is computed to be 6.5 /iK, is subtracted from this in quadrature 
then one finds an "intrinsic" error of 8.2 /iK. The full data set should have an intrinsic 



error of 8.2A/l30h/428h fiK = 4.5 /iK which is close to the computed value of 5.5 /iK. In 



^The slight differences between Peebles (1994) and Bond (1994) may be ignored for large 
multipoles, l. 



addition, by inspection of Figure 2 and Table 1, one sees that the data from the three 
radiometers may differ by more than one would expect from random noise alone, possibly 
due to the atmosphere. This would tend broaden L{Arms, P)- 

There are two sources of a systematic error on /?, uncertainty in the relative calibration 
between channels and the uncertainty in the fit to the Cas-A calibration data . These 
combine to produce an error of ±0.08. In addition, we cannot be certain that atmospheric 
fiuctuations are not biasing our answer. While the bottom two entries in Table 1 agree 
statistically, one might worry that low cut levels result in more negative indices. To ensure 
that the quoted value for j3 represents one's confidence, we quote the average value of f3 
with an error bar that nearly encompasses the errors from both entries. Including the 
systematic uncertainty, P = — O.llgJ. 

To investigate the possibility that Arms is due to a Galactic foreground, we examine 
L{Arms, P) bearing in mind that there might be an atmospheric bias. L{Arms, P) is 
normalized to 1 at its maximum. For free-free emission, max[L{Arms, ~'^-^)] = 2 x 10~^ 
and 3 x 10~^ for the bottom two entries in Table 1 respectively. For interstellar dust 
emission, maxlLlArms, 1-6)] = 4 x 10"^ and 2 x 10~^ for the same two entries. Clearly 
these possibilities are strongly ruled out, in the worst case at 3o" {L = 0.003 for a normal 
distribution). In Wl, Arms from extra-galactic sources is expected to be ^ 9 /iK. A 
conservative index for these sources is /9 ~ —2; thus, this explanation too is disfavored. An 
anisotropy in the CMB, /9 = 0, is the most compelling interpretation of the data. 

In the analysis, we assume that only one celestial component produces the anisotropy. 
Dodelson & Kosowsky (1994b) have performed a more general analysis on the SK93 data 
that allows two fiuctuating components, CMB and free-free. The result was a decrease in 
the significance of the detection of CMB. With this new limit on (3, Arms of uncorrelated 
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free-free emission [| is less than ^ 27 /iK, an additional constraint on two component models. 
Using these data alone, we cannot rule out a model which assumes co-located interstellar 
dust and free-free emission, but, based on FIRS (Ganga 1994) we expect dust to contribute 
less than 1 /iK, again, added in quadrature. 

Gaussian correlation functions are often used to compare measurements. We find 
Co = 60tJ^ for a correlation angle 9c = 1.0°. Alternatively, the result may be interpreted 
in the context of the 'standard' CDM model (Steinhardt 1993; Bond et al. 1994), 
{Qf, = 0.05, Hq = 50 km/s/Mpc,A = 0, n^ = 1, fi = 1, no reionization), for which we find 
that the root-mean-squared quadrupole has a value Q = 18.9l3g /^K. 

We would like to thank George Sofko and Mike McKibben at the The University of 
Saskatchewan and Larry Snodgrass at the Canadian SRC for their valuable assistance 
during the 1994 observing season. We are indebted to Marian Pospieszalski and Mike 
Ballister of NRAO who provided the amazing HEMT amplifiers that made this experiment 
possible. We also thank Chris Barnes, Carrie Brown, Weihsueh Chiu, Randi Cohen, Peter 
Csatorday, Cathy Cukras, Peter Kalmus, John Kulvicki, Wendy Lane, Young Lee, Naser 
Quershi, and Peter Wolanin for building many of the components of the experiment. This 
work was supported by NSF grant PH 89-21378, NASA grants NAGW-2801 & NAGW-1482, 
and and NSF NYI grant to L. Page. 

The data, analysis software, beam profile, and correlation coefficients will be made 
publicly available upon acceptance of this Letter. 



^The free-free (or dust) contribution is estimated by assuming that it adds in quadrature 
to the CMB signal. We assume acMB/crff =if3ff—f3)/il3cMB—(3). Then, A^^^ = crcMB+^}f- 
For Arms = 44 /iK and [3 = —0.9, ajf = 27 /iK, and (Tcmb = 35 /iK, near the bottom of our 
confidence range. 
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Table 1. Summary of Results 



Data Set Total time^ C'' Mean offset" A^„/ (3 5Tf 



(h) (mK/deg) (/iK) (/iK) (^K) 



Ka93f 130 4.5 -60/ -360 

79 3.0 

Ka93, A-Bs. 130 4.5 

Ka94 140 2.5 +110/ -340 

107 1.7 

Ka94, A-Bs 140 2.5 

Q94 158 2.5 ■ ■ ■ / - 400 

98 1.7 

Ka93-Ka94 298 4.5,2.5 

Ka93 + Ka94 298 4.5, 2.5 

186 3.0,1.7 

Ka93 + Ka94, A-Bs 298 4.5, 2.5 

Ka94 + Q94 298 2.5,2.5 

205 1.7,1.7 
Ka93 + Ka94 + Q94 428 4.5, 2.5, 2.5 

284 3.0, 1.7, 1.7 



371^2 


-0.28l?ii 


361^^ 


44+i^ 


-o.69t?:I^ 


42+lg 


0,31 






49+^2 


-0.701°;?^ 


AltY 


521^^ 


-o.77i°:i 


52tf 


0,41 






6lli^ 


+i.35t^:^4 


57tl? 


52+il 


-5.97f,f 


4811? 


0,33 






44+^^ 


-0.43+^:^^ 


A2tf 


49i^2 


-0.62l°;i 


AltY 


0,34 






48+^^ 


+o.ioi°:^^ 


AAtf 


511^^ 


-o.54t°:i 


Alt]' 


AAtf 


+0.17tn^ 


A2tf 


Alt]' 


-0.40+°:^^ 


AAtf 
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" Between February 21 and March 15 of 1993, the sky was observed with the Ka93 
radiometer for 450 hours. Between January 13, 1994 and February 28, 1994, the sky 
was observed for 243 hours with the Ka94 radiometer data, and 363 hours with the 
Q94 radiometer. 

^ The cut level, ( is the root-mean-squared amplitude in the horizontal temperature 
gradient (single-difference) for Is of data. The larger cut levels correspond to factors 
of 3.1, 2.5, and 2.5 above system noise (Tgys = Tree + Tatm) noise for the Ka93, 
Ka94, and Q94 radiometers respectively. The single difference in the 1994 observing 
scheme is much more sensitive to atmospheric gradients than that from 1993. Data 
from the single-difference are not used in the anisotropy analysis. The data from 
the double-difference, which we analyze for anisotropy, are dominated by the system 
noise. 

'^ The offsets are given as (vertical polarization. A) /(horizontal polarization, B). 

'^ Data are quoted as the maximum of the likelihood with the range of the error bars 
encompassing the 16% to 84% integral of the likelihood marginalized with a uniform 
prior. Where there is no "2(t" lower bound, the maximum and 95% upper bound 
are given. The units are thermodynamic temperatures relative to a 2.726 K Planck 
emitter. At 30 GHz, dTant/dT = 0.98; at 45 GHz, dTant/dT = 0.95. Given our finite 
sampling of the CMB, presumed to be a random field, it is not possible to determine 
Arms to better than ^ ± 6.5 fiK. 



The synthesized beams differ by about 4% between channels. For Ka93 the window 

eff 



spans leff = 66^20 ^^(^ \/-^(^) = 1-04 (Bond, 1994); for Ka94 the window spans 



/^^^ = 67^20 and Jl{W) = 1.04; and for Q94 the window spans /e// = 71^1 and 



I{W) = 1.08. Arms = STi/Jl{W). For the Ka and Q combination, we use 
leff = 69^22 and Jl{W) = 1.06. 
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^ The values for the Ka93 data have been recomputed given our increased understanding 
of Cas-A and its environs. The original values were A^ms = 33 ± 10 /iK and 
P = — 0.3ti;2- The cut level ( = 4.5 mK/deg corresponds to r] = 3 in Wl. 

^ This is a measure of the anisotropy of the difference between two perpendicular 
polarizations. It is also a gauge of the consistency between the A and B channels. 
The values are for (3 = 0. 
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13^ 12^ ^^^ ^Q^ 



Fig. 1. — The region of observation and beam pattern for SK93. The contours show the 
beam profile for Is of data. In the analysis, all the data with the beam center within one 
hour of right ascension are combined; thus, the above beam pattern is smeared out. Note 
how much the pattern resembles a "double-difference" (positive lobe: solid lines; negative 
lobes: dashed lines). The '+' signs mark the sources in Kiihr et al. (1981). Stars mark 
the sources with flat spectra. The size of the symbol is proportional to the flux. Lines of 
constant galactic latitude are also shown. Note the slight curvature of the observing pattern 
with respect to lines of constant longitude; this is a result of the azimuthal chop. Note that 
if there were a hot spot in the CMB near 5 = 87° and a = 5^, the data in bins 4 through 8 
would tend to be low. 
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Fig. 2. — The data for Ka93, Ka94, and Q94 as a function of right ascension bin. Data from 
the six Ka bands have been added together as have the data from the three Q bands. Also, 
the data taken in the east have been combined with the data taken in the west. The error 
bars include the effects of the correlations between channels. There are no correlations, other 
than those from the CMB, between the three data sets. Because of the observing strategy, 
each bin is correlated with the adjacent bins. All correlations have been accounted for in 
the analysis. 



